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The largest Lyapunov exponent has been examined in the dynamical-mode locking phenomena of 
the ac-|-dc driven dissipative Frenkel-Kontorova model with deformable substrate potential. Due to 
deformation, large fractional and higher order subharmonic steps appear in the response function 
of the system. Computation of the largest Lyapunov exponent as a way to verify their presence led 
to the observation of the Farey sequence. In the standard regime, the appearance of halfinteger and 
other subharmonic steps between the large harmonic steps, and their relative sizes follow the Farey 
construction. In the nonstandard regime, however, the halfinteger steps are larger than harmonic 
ones, and Farey construction is only present in the appearance of higher order subharmonic steps. 
The examination of Lyapunov exponents has also shown that regardless of the substrate potential 
or deformation, there was no chaos in the system. 

PACS numbers: 05.45.-a; 45.05.-|-x; 71.45.Lr; 74.81.Fa 


I. INTRODUCTION 

In the examination of Shapiro steps, finding the best 
method or tool to verify their presence has been the mat¬ 
ter of many studies in various physical systems. Numer¬ 
ous theoretical and experimental results on Shapiro steps 
obtained in dissipative systems such as charge- or spin- 
density waves conductors the systems of Joseph- 

son function arrays and superconducting nanowires 
[^, 0 have been the main impulse and motivation for 
our studies of the ac-|-dc driven overdamped (dissipative) 
Frenkel-Kontorova (FK) model 0- It is well known that 
when these systems are subjected under an external ac 
driver, their dynamics is characterized by the appearance 
of Shapiro steps. These steps are due to interference phe¬ 
nomena or dynamical mode-locking (synchronization) of 
the internal frequency with the applied external one. If 
the locking appears at the integer values of the exter¬ 
nal frequency, the steps are called harmonic while for the 
locking at rational (noninteger) values of frequency they 
are called subharmonic. 

In realistic systems due to presence of noise, impuri¬ 
ties and other environmental effects, detection of Shapiro 
steps, particularly the subharmonic ones, is usually very 
difficult. On the other hand, in theoretical works, their 
observation could also be a problem since their size is 
often so small that they are invisible on the regular plot 
of the response function. One of the most sensitive ways 
to verify the existence of Shapiro steps is the calculation 
of the largest Lyapunov exponent [ll|. Always when the 
system is dynamically mode-locked, the largest Lyapunov 
exponent has negative values [l2, [0 ■ Therefore, an ex¬ 
amination of the largest Lyapunov exponent for some 
interval of driven force will precisely reveal the presence 
of any harmonic or subharmonic mode-locking. 


Calculation of the largest Lyapunov exponent has been 
already used as a way to examine the existence of subhar¬ 
monic Shapiro steps in the standard FK model 00. 
The standard Frenkel-Kontorova (FK) model represents 
a chain of harmonically interactin g p articles subjected 
to a sinusoidal substrate potential 0. It describes dif¬ 
ferent commensurate or incommensurate structures that 
under an external driving force, show very rich dynami¬ 
cal behavior. In the presence of an external ac-|-dc driv¬ 
ing force, the dynamics is characterized by the appear¬ 
ance of the staircase macroscopic response or the Shapiro 
steps in the response function v{F) of the system 0-[0. 
Though the standard FK model has been very successful 
in the studies of some effects related to Shapiro steps, 
its applications is still very restricted. Namely, in the 
standard FK model, the subharmonic steps either do not 
exists in the case of commensurate structures with in¬ 
teger values of winding number 00 or their size is 
so small that analysis of their properties is very difficult 
0 - [T3| . The absence of subharmonic steps for the com¬ 
mensurate structures with integer value of winding num¬ 
ber, and their presence in the case of rational (noninte¬ 
ger) winding number was confirmed by the calculation of 
the largest Lyapunov exponent [0 . However, contrary 
to the standard case, the large subharmonic steps can 
appear in any commensurate structure of the nonstan¬ 
dard FK model such as the one with the asymmetric de¬ 
formable substrate potential (ASDP) 0. This potential 
belongs to the family of nonlinear periodic deformable 
potentials, introduced by Remoissent and Peyrard 0 
as the way to model many specific physical situations 
without employing perturbation methods. 

In this paper, by using the largest Lyapunov expo¬ 
nent computation technique, we will examine the ap¬ 
pearance of both harmonic and subharmonic steps in the 
FK model with asymmetric deformable substrate poten- 
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tial (ASDP). In the analysis of the largest Lyapunov ex¬ 
ponent, we have observed one interesting property: the 
Shapiro steps and their relative sizes appear according 
to the Farey construction only in the standard regime 
when large harmonic steps are dominant in the response 
function. The paper is organizes as follows. The model 
and methods are introduced in Sec. II, and the results 
are discussed in Sec. III. Finally, Sec. IV concludes the 
paper. 


II. MODEL AND METHODS 

We consider the dissipative (overdamped) dynamics of 
a series of coupled harmonic oscillators ui subjected in 
the ASDP: 

V(u) = ^ (l-r^)Ml-cos(27ru)] 

(27r)2 [i + ^2 + 2rcos(^u)]^ ’ 

where K is the pinning strength, and r is the shape or 
deformation parameter (—1 < r < 1) which can be var¬ 
ied continuously. By changing the shape parameter r the 
potential can be tuned in a very fine way, from the simple 
sinusoidal one for r = 0 to deformable for 0 < r < 1. In 
Fig[Tl the commensurate structure w = 1/2 in ASDP is 




FIG. 1: (Color online) Particles moving in asymmetric de¬ 
formable potential for w = |, A = 4, and two different values 
of the shape parameters r = 0 and r = 0.5. Particles are 
represented by red dots. 

average interparticle distance lu = — ui), or the 

so called winding number is one of the main parameters 
that describes the FK model. The system exhibits com¬ 
mensurate phase for rational values of winding number 
w, and incommensurate phase for irrational ones. 


In the present paper, the system of coupled harmonic 
oscillators in ASDP is driven by the dc-|-ac forces F{t) = 
F + FacCOs{2TTi/Qt). The equations of motion are 


ill = ui+i + ui-i - 2ui 


dV 


( 2 ) 


where I = —-j; Y’ position of Zth particle, F 

is dc force, where Fac and 27rr'o are the amplitude and 
the circular frequency of the ac force, respectively. Since 
the substrate potential is homogeneous (it does not de¬ 
pend of the particles index 1) relabeling of the position 
of particles will not change the properties of the config¬ 
uration 0- When the system is driven by an external 
ac-|-dc force, two different frequency scales appear in the 
system: the frequency of the external periodic force vq, 
and the characteristic frequency of the motion of parti¬ 
cles over the ASDP driven by the average force F. The 
competition between those frequency scales can result in 
the appearance of resonance (dynamical mode locking or 
Shapiro steps). 

Solution of the system ([2]) is called resonant if average 
velocity v satishes the relation: 


ioj + ja 

V = -l/Q, 

m 


( 3 ) 


where triplet (i^j,m) are integers numbers. Resonant 
velocity is called harmonic if to = 1 and subharmonic if 
TO 7 ^ 1. (In case of w = | we can use v = for 

marking harmonic and subharmonic steps). Parameter a 
is the period of the potential V(u) and in the case of no 
deformation a = 1, and with deformation a = 2 as can 
be seen in Fig[T] For a rational value of w = p/q {p and q 
coprime integers) the triplet is not unique (this triplet is 
unique only for incommensurate structures 0 ). In this 
paper we will consider only the commensurate structure 
w = 1/2. 

The equations of motion ([2]) have been integrated by 
using a fourth order Runge-Kutta method with periodic 
boundary conditions. The time step used in simulations 
was O.Olr, where t = —. The force is varied adiabati- 
cally with the step 10“®. 

We shall be focused on calculating the largest Lya¬ 
punov exponent A [l^ . It is well known that the Lya¬ 
punov exponent gives a quantitative measure on the pres¬ 
ence of chaos in dynamical systems 0 , however, it also 
proves to be extremely sensitive to the existence of both 
harmonic and subharmonic Shapiro steps. When the sys¬ 
tem is dynamically mode locked, i.e. on the step, the 
trajectories of particles are periodic in time which is re¬ 
flected by the negative value of the largest Lyapunov ex¬ 
ponent. Outside the steps, where there is no onset of dy¬ 
namical mode locking, the trajectories are quasi-periodic 
which is confirmed by the zero of the Lyapunov exponent 
( 00 ). We choose an appropriate perturbed point u[ 
in our computations according to: 


u'litss) 


Ui{tss) ± 



( 4 ) 
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where tss is time when the steady-state has been achieved 
in our system, and dp is the parameter that expresses the 
change in the initial positions of particles of the model. 
In order to make sure that projecting is always done 
onto the subspace dominated by the largest Lyapunov 
exponent, the sign in front of the square term in Eql4] 
is randomly selected where the plus and minus sign ap¬ 
pear with the same probability. We sample and readjust 
following Sprott every 25 or so time steps. In our 
calculations, we used tss = 300 t and do = 10“”^. For 
convenience, in further text, the largest Lyapunov expo¬ 
nent will be denoted just as the Lyapunov exponent. 


III. RESULTS 

In the present paper, the Lyapunov exponent is exam¬ 
ined for different deformations of the substrate potential. 
In Fig. [21 the Lyapunov exponent as a function of the 
driving force for three different values of deformation pa¬ 
rameter is presented. The insets show the correspond¬ 
ing response functions v{F) (the average velocity as a 
function of average driving force). As one can see, the 
computed Lyapunov exponents are always A < 0, which 
implies that with the change of deformation r we have 
not introduced chaos in our system (presence of chaos 
would result in the positive values of Lyapunov expo¬ 
nent). Domain of F in Fig. |21 for which we calculated 
the exponent, differs with r due to the fact that for differ¬ 
ent values of r the same steps appear in different regions 
of F (see [I3,[l^)- In the standard case in Fig. 121(a), we 
can see the large minima which correspond to harmonic 
steps and which size changes monotonically. As defor¬ 
mation increases in Fig. 121(b) and (c), the minima which 
corresponds to the large halfinteger and higher order sub¬ 
harmonic steps appear where the changing of their size 
is not monotonic any more. 

Using the Fq.®, the Shapiro steps could be now iden¬ 
tified. It is well known that in the standard FK model 
(r = 0) with integer value of winding number, there 
would be no subharmonic mode locking 0 , and con¬ 
sequently, no steps between harmonic ones on the plot 
of response function v{F). On the other hand, when 
w = 1/2, only halfinteger step 3/2 which appears be¬ 
tween the first and the second harmonic could be visible 
[T^ . However, computation of the Lyapunov exponent 
between first and the second harmonic steps reveal other 
subharmonic steps as can be seen in Fig. [31 The areas 
under the minima correspond to the size of the steps, 
i.e. for larger step, the area under the minimum will 
be larger. If we examine the subharmonic steps in Fig. 
121 we could see that the first largest fractional step be¬ 
tween the step I and the step 2 is the step 3/2. Then, 
the largest step between the steps I and 3/2 would be 
the step 4/3 while the largest one between the steps 3/2 
and 2 is 5/3. Therefore, according to the appearance of 
fractional steps between the first I/I and the second 2/1 
harmonics we may write the following sequence: 



F 



F 



F 


FIG. 2: (Color online) The Lyapunov exponent as a function 
of the average force for commensurate structure oj = 1/2, K = 
4, Fac = 0.2, Vo = 0.2 and three different values of shape pa¬ 
rameter (a) r = 0, (b) r — 0.28 and (c) r = 0.6. The insets 
show the corresponding response functions v{F). 


1 4 3 5 2 

T’3’2’3’1' 


(5) 


This sequence of numbers represents the Farey se¬ 
quence well known in number theory pm. 

The Farey sequence Fjq- of order J\f is an ascending 
sequence of irreducible fractions between 0 and 1, whose 
denominators are less or equal then Af (mill. The first 
few would be: 


Fi 

Tz 

To 


/ 0 

ll 









11’ 

IJ 









r 0 

1 


1 







1 ^ ’ 

2 ’ 

1 








1 0 

1 

1 

2 


1 





1 1 ’ 

3’ 

2 ’ 

3’ 

1 

r 





1 0 

1 

1 

1 

2 

3 


1 



1 1 ’ 

4’ 

3’ 

2 ’ 

3’ 

4’ 

1 

r 



( 0 

1 

1 

1 

2 

1 

3 

2 

3 

4 

1 ^ ’ 

5’ 

4’ 

3> 

5> 

2’ 

5> 

3’ 

4’ 

5’ 1 / 


( 6 ) 








































4 



F 

FIG. 3: (Color online) Lyapunov exponent as a function of 
average driving force for r = 0 (the rest of parameters are as 
in Fig. II). The inset shows the response function v{F) drawn 
for the same interval of force. This result is obtained in [l^. 


Therefore, if we have two rational fractions in Farey se¬ 
quence - (jj, q are coprime integers) and ^ ip' ,q' are co¬ 
prime integers), the rational fraction which lies between 
them and has the smallest denominator is 


q” q q' q + q' 


(7) 


where p", q" are coprime integers. This statement is 
trivially extended to the case of interval between 1 and 2 , 
and further on (Theorems 28 and 29 in [l^). The largest 
step between | and ^ 7 , if exists, will be step | © ^ 7 . 
The Farey sequence could be easily understood from the 
diagram in Fig. |4] 



FIG. 4: (Color online) Section of the Farey construction (rep¬ 
resented as a rooted binary tree graph) a) from 0 to 1, b) from 
1 to 2 according to Eq. [3 


For example, in the case of the FK model with the in¬ 
teger value of winding number, there is no subharmonic 
mode locking, which implies there is only the Farey se¬ 
quence of order one. However, if the winding number 
is rational noninteger such as the case w = 1/2 in Fig. 
[21 one can see that the largest step between harmonic 


steps 1 and 2 is halfinteger step 3/2. From the set the¬ 
ory ( 2 ^ we know that between any two rational fractions 
lie countable many, Hq rational fractions and therefore, 
countable many possible Shapiro steps between any two 
harmonic steps in our model. 

If the potential gets deformed, the large halfinteger 
step and higher order subharmonic steps appear [17], [l^ . 
Contrary to the case r = 0 in Fig. |3l now for r = 0.01 in 
Fig. O the large 4/3 and 5/3 steps are clearly visible. The 



FIG. 5: (Color online) Lyapunov exponent as a function of 
average driving force for r — 0.01 (the rest of parameters are 
as in Fig. 13). The inset shows the response function v{F) 
drawn for the same interval of force. 

higher order subharmonic steps, such as 4/3 and 5/3 (to 
the left and to the right), are appearing in a symmetric 
manner with respect to the step 3/2. 

With the further increase of deformation r, the step 
widths increase faster on the right side from halfinteger 
step 3/2 than on the left one as can be seen in Fig. |6l 



FIG. 6: (Color online) Lyapunov exponent as a function of 
average driving force for r — 0.05 (the rest of parameters are 
as in Fig. m. The inset shows the response function v{F) 
drawn for the same interval of force. 

It was shown in our previous work , that size of 
halfinteger and subharmonic Shapiro steps increase with 
deformation and after reaching their maxima for some 
value of r, decrease to zero. If we calculate the Lyapunov 
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exponent for r = 0.28, which is the value of r for which 
halfinteger step reaches its maximum, we obtain results 
presented in Fig. 0 At this value of deformation some 



FIG. 7: (Color online) Lyapunov exponent as a function of 
average driving force for r = 0.28 (the rest of parameters are 
as in Fig. [2|. The inset shows the response function v{F) 
drawn for the same interval of force. 

higher subharmonic steps already start to disappear. 

At large deformation of the potential, the size of halfin¬ 
teger steps decreases, and the most of higher order sub¬ 
harmonic steps have completely vanished [T^ . This is 
conhrmed by the results in Fig. |51 where the Lyapunov 
exponent for r = 0.5 has been calculated. Disappearance 



FIG. 8: (Color online) Lyapunov exponent as a function of 
average driving force for r = 0.5 (the rest of parameters are 
as in Fig. [2]). The inset shows the response function v{F) 
drawn for the same interval of force. 

of steps is also clearly visible in Fig. [2j 

Deformation of the potential obviously strongly affects 
the steps as we can see in Fig. IBUS] It appears that with 
the increase of the deformation r, the right side of the 
Farey construction is heavily favored over the left one. 
In particular, we observe that at each level of our binary 
tree graph (Farey construction in Fig. 2]) the child node 
(or step) that takes preference is the one on the right. 
This means that with the increase of the deformation the 
steps that are present and become increasingly dominant 
are 3/2,5/3, 7^, 9/5. 

We have analyzed also the systems with other types of 


deformable potentials [2^, such as variable, double bar¬ 
rier and double well potential, and we have been always 
able to observe the appearance of steps in accordance 
with the Farey construction [HI- Therefore, for two steps 
I and ^, the next largest step between them will be 
where denominator determines the size of steps in terms 
that the size of steps decreases as the denominator in¬ 
creases. It is important to note that Farey construction 
tell us the order and the relative sizes of steps but it does 
not tell us the actual step width or why they appear [Hi¬ 
lt is well known that the size of harmonic and halfin¬ 
teger steps are correlated, whereby the larger the size 
of harmonic the smaller the one of halfinteger step and 
vice versa In some cases, depending on the 

system parameters, the size of halfinteger steps could be 
even larger than the size of harmonic ones jl9| . The size 
of halfinteger and other subharmonic steps strongly af¬ 
fects the behavior of the system, and according to that in 
previous works [l^ [2^ , the three different types of sys¬ 
tem behavior have been classified: the standard behavior 
for small halfinteger steps, the behavior for intermediate 
halfinteger steps and the behavior in the presence of large 
halfinteger steps. 

If we have two harmonic steps, according to Farey se¬ 
quence the next largest step which appears between them 
is the halfinteger stem but this is not the case for nonstan¬ 
dard behavior [l^ , since halfinteger steps are larger 

than harmonic ones. In such case, could we still have 
the presence of Farey sequence? In Fig. [31 the response 
function in the case of large halfinteger steps is presented. 
It is obvious that in the nonstandard case, the relative 



FIG. 9: The average velocity as a function of average driving 
force for Fac = 0.55, and r — 0.2 (the rest of parameters are 
as in Fig. o. The numbers mark halHnteger and subharmonic 
steps. 

sizes of harmonic and halfinteger steps do not follow the 
Farey construction, and going from harmonic to halfinte¬ 
ger steps, the size of step does not decreases as denomi¬ 
nator increases, on the contrary, the halfinteger steps 1/2 
and 3/2 are larger than harmonic ones 1/1 and 2/1. How¬ 
ever, the higher order subharmonic steps between halfin¬ 
teger and harmonic steps still appear according to Farey 
construction and their sizes decrease as the denominator 
increases. 
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Calculation of the Lyapunov exponent gives a possibil¬ 
ity not only to detect all resonances in the response func¬ 
tion, but also to detect the presence of chaos. In all our 
simulations performed on the ac-|-dc driven overdamped 
FK model we did not observed any chaos. Presence of 
deformable substrate potentials and different level of de¬ 
formations did not introduced chaotic behavior into the 
system. Contrary to our case, chaos has been observed 
in the spatiotemporal dynamics of moving kinks in the 
damped dc driven FK model where the resonances ap¬ 
pear due to competitions between the moving kinks and 
their radiated phase modes [2^. Also, structured chaos 
has been observed in a Josephson junction systems where 
chaotic regions appear between the subharmonic Shapiro 
steps at certain values of system parameters [^ . 

IV. CONCLUSION 

In this paper, we have presented detailed analysis of 
the Shapiro steps in the ac-|-dc driven dissipative FK 
model by using the Lyapunov computation technique. 
The obtained results show the presence of Farey sequence 
in the appearance of subharmonic steps. The steps and 
their relative sizes follow exactly the Farey construction 
only in the standard regime when harmonic steps are 
the largest one. However, in the nonstandard regime, 
the halfinteger steps are larger than harmonic ones, and 
Farey sequence appears only in the order and relative 
sizes of higher order subharmonic Shapiro steps. Lya¬ 
punov exponent analysis is certainly one of the best ways 
to get an accurate answer about the presence of chaos in 


the system. Computations of Lyapunov exponent have 
been performed for different system parameters, and re¬ 
gardless of the deformation, no chaos has been observed 
in the behavior of the system. Absence of chaos in the 
presence of deformable potentials certainly requires fur¬ 
ther investigation. This problem and the possibility of 
chaotic behavior in other situations such as presence of 
noise will be subject of future examinations. 

Presented results could be important for the studies of 
Shapiro steps in all systems that are closely related to the 
dissipative dynamics of the FK model. In experimental 
and theoretical works performed in charge density wave 
systems and the systems of Josephson junction arrays, 
measuring of differential resistance is usually used to de¬ 
tect subharmonic steps. If we look for example at the 
results obtained in sliding charge-density wave systems 
0, the systems of Josephson junction arrays @,0 or su¬ 
perconduction nanowires Q , we can observe the presence 
of Farey construction in the appearance of Shapiro steps. 
Our analysis shows that Farey construction can not be 
always generally applied when it comes to relative sizes 
of the observed steps. Since the appearance and origin of 
the subharmonic Shapiro steps are still a matter of de¬ 
bates, we hope that these results could bring an insight 
into understanding of these physical phenomena. 
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